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Abstract. Although Casimir, or quantum vacuum, forces between distinct bodies, or 
self-stresses of individual bodies, have been calculated by a variety of different methods 
since 1948, they have always been plagued by divergences. Some of these divergences 
are associated with the volume, and so may be more or less unambiguously removed, 
while other divergences are associated with the surface. The interpretation of these has 
CO , been quite controversial. Particularly mysterious is the contradiction between finite 

total self-energies and surface divergences in the local energy density. In this paper we 
clarify the role of surface divergences. 
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1. Introduction 



The subject of local energy density associated with the confinement of quantum fields 
by surfaces has a rather long history. For example, Brown and Maclay pQ computed 
the vacuum expectation value of the electromagnetic energy-momentum tensor between 
two parallel perfectly conducting plates, which is twice that of a conformally coupled 
massless scalar field satisfying Dirichlet or Neumann boundary conditions on the plates, 
namely for plates separated by a distance a in the x direction, 

(^ = 77^^ (-1,-3, 1,1), (1) 

which corresponds precisely to the attractive energy or pressure found by Casimir [2] in 
the same situation. If a nonconformal scalar stress tensor is used, a position-dependent 
term in the stress tensor appears, which does not contribute to either the total energy 
or the pressure on the plates [HI Ej . 

Local surface divergences were first discussed for arbitrary smooth boundaries by 
Deutsch and Candelas J5j. They found cubic divergences in the energy density as 
one approaches the surface; for example, outside a Dirichlet sphere (that is, for a 
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Figure 1. A plot of the TM Casimir stress F™ for — 2 < D < 4 on a spherical shell, 
compared with F TE , taken from Bender and Milton |8] For D < 2 (D < 0) the 
stress F™ (F TE ) is complex and we have plotted KeF. 



conformally-coupled scalar field satisfying Dirichlet boundary conditions on the surface) 
the energy density diverges as 

1 1 

r — > a+ : u ~ tt—, r=-, 2 

360vr 2 a(r-a) 3 ' V ; 

where a is the radius of the sphere. 

This raises the question: How can it be that the total Casimir energy of a 
Dirichlet sphere (or a perfectly conducting sphere in electrodynamics) is finite? The 
electromagnetic case is the well-known one first calculated by Boyer jH], £ EM = 
0.04618/a, while the scalar case was first worked out by Bender and Milton [Jj, 
S s = 0.002817/a. In general the Casimir energy of a region bounded by a perfect 
hyperspherical surface depends in a complicated way upon the number of spatial 
dimension D, as shown in Figure d 

Thus there has been a suspicion since the time of Deutsch and Candelas that 
there was something incomplete in the calculations of Casimir self energies of ideal 
closed boundaries. (We note that there is now a proof that any such smooth perfectly 
conducting boundary possesses a finite electromagnetic Casimir energy Whether 
such an idealized limit is physical is, of course, another question.) This suspicion has 
been recently intensified by a series of talks and papers by Graham et al ^Uj- The 
essential outcome of their analysis is that for a (^-function sphere, described by the 
following Lagrangian for a massless scalar field, 



2 M 



IA 

2a? 



5(r — a) 



(3) 



a divergence occurs in third order in A. (They claimed a divergence in second order, 
but that was spurious, arising from the omission of a surface term in the integration 
by parts [SI HU [T2] See (J3U|) below.) This divergence in 0(A 3 ) in fact was discovered 
much earlier by Bordag, Kirsten, and Vassilevich [T5| , and possible ways of dealing with 
it have been suggested [TJJ[T5]. Objections complementary to those of Graham et al 
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have also been voiced by Barton [TBI AT] , all of which raise doubts as to the physical 
relevance of results such as those of Boyer. 



2. Green's function for A sphere 

We consider the potential 
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hit 



2a 2 



2 cr(r), where cr(r) 



r < a_, 

a_ < r < a + , 

a + < r. 



(4) 



Here a± = a ± 5/2, and we set ft5 = 1. We have chosen the dimensions of A so that 
the total energy of interaction does not explicitly refer to the radius a. In the limit as 
5^0 (or h — > oo) we recover the 5-function sphere. 
The Green's function equation, with k 2 = —u> 2 , 



k 2 + ^o] £(r,r') = 5(r-r') 



(5) 



may be straightforwardly solved. We introduce the reduced Green's function 

£(r,r') = Y,9i{r,r')Y lm {e,<l>)YL{0'A'), (6) 

Im 

which in turn may be expressed in terms of the modified Riccati-Bessel functions, 



si{x) 




K l+1/2 (x). 



(7) 



The reduced Green's function is, outside of the shell: 



1 



r, r < a_ 



r, r > a + : 
Here the denominator is 



9l 



9l 



1 

nrr' 



sj(ftr<)e/(/«r>) - ^ei(Kr)ei(Kr') 



(8b) 



H = [KS;(Ka_)e;(fi;'a 



- K's;(Ka_)e;(fi;'a_)] [/t'e;(fi:a_|_)s;(fi; / a + ) — /cej(/ca + )sj(/c'a+)] 

— [/tS;(/ta_)si(/t'a_) — K,'si(K,a-)s'i(K,'a-)] [K / e;(Ka + )eJ(/t / a + ) — /tej(/ta + )e; (fi/a+)], (9) 

while the numerator H is obtained from H by replacing sj(/«a_) — > e/(/ta_), and 5 
is obtained from H by replacing e^(/«a + ) — > sz(/ta + ). Here = V^ 2 + Aft. Green's 
function within the shell, a_ < r < a + , is given by 



5/ = -; — ri si(K'r < )ei(K'r > ) - — [s i (« / r)e J (/cV') + si(K?r')ei(K'r)] 
Krr I ^ L 

x [KeJ(/«a + )ei(fi; / a + ) — K'ez(Ka + )e|(K^i+)][ttS;(Ka_)sz(fi/a_) — K / s/(/«a_)sJ(«; / a_)] 

— s;(/tV / )si(K / r)[KeJ(Ka + )ei(fi; / a + ) — /t / e/(Ka + )eJ(K / a + )] 
x [KsJ(/ta_)e;(fi: / a_) — K / s/(Ka_)eJ(/« / a_)] 

— ei(K ! r')ei(K ! r) [/€eJ(Ka+)s/(K / a + ) — /t'e/(/ta + )sJ(K / a + )] 

x [KsJ(/ta__)s/(fi: / a_) — K'si(Ka-)s'i(K'a-)] 



(10) 
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3. Energy density 

We can calculate the local energy density from the stress tensor: 

= d^d^ - g^C - £(d>*d" - g^d 2 )^ 2 , 
from which the energy density follows: 



& 



2 J.2 



(12) 



where the conformal value is given by £ = 1/6. To obtain the one-loop vacuum 
expectation values, we use the connection to the Green's function 



The energy density thus is, within or outside the shell, 
(T 00 ) = i ( d°d'° + V ■ V + { Xh/a2 )) G{x, x') 



iv 2 G(x,x). 



(13) 



(14) 



When we insert the partial wave decomposition of the Green's function (JBjl. the 
expression for the energy density is immediately reduced to (inside or outside the shell, 
but not within it) 



(T o 0) = r 21 + 1 

[ ' Jo 2tt^ 4tt 



-k + d r d r i + 



. 1 <9 2 <9 . . 
2£^^-r — ^(r,r) 



(15) 



We insert the Green's function in the exterior region, but delete the free part, the 
first term in (|8ojl . (j8fej) . (jlUjl . which corresponds to the bulk energy which would be 
present if either medium filled all of space, leaving us with for r > a + (for r < a_, 
ei -> s/ and S ->• H) 



n(r) = -(l-40/ 



2Z + 1H fe? 



47T 5 



.KT) 



2 l + 4£ /(/ + !) 1 



4£ 



- e/ (Kr) e ;(/tr) + — ef («r) k 



(16) 



4. Surface divergences 

We want to examine the singularity structure as r — > a + +. For this purpose we use the 
leading uniform asymptotic expansion, I — > oo, 



e 



1 1 



zt 



Sl {X) 



2Jzi 



e»\ (17) 



where u = 1 + 1/2, a; = vz, t = (1 + z 2 )~ 1/2 , drj/dz = 1/zt. 

Let us consider the thin shell limit, 8 — *■ 0, h5 — 1, where it is easy to check that 
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which is exactly the coefficient occurring in the (^-function potential. There are two 
simple limits of this, strong and weak coupling: (na ~ 1) 

-^oo: ---Li, _-0: --_ s ?( Ka . (19) 
a n ei[K,a) a s «a z 

In either case, we carry out the asymptotic sum over angular momentum using the 

uniform asympotic expansion and 



^ 2 sinh % ' 

i=o z blml 2 



r](z) — rj (z— 



,, r — a 2 r — a , , 

2«/(z)— = 7 — . (20) 



-^oo: it T7 — -, (21a) 



n) '■ " ^ 4_n 



0: U W~ (--) ^7(1-60(^1 , n<4, (216) 



The remaining integrals over z are elementary, and in this way we find that the leading 
divergences are as r — > a+, 

A 1 1 — 6^ 

> OO : It ~ -7 r 

a 16n 2 (r — a) 

a \ a J 96tc a 4 

the latter being the leading divergence in order n, which clearly seems to demonstrate 

the virtue of the conformal value of £ = 1/6. (The value for the Dirichlet sphere first 

appeared in Deutsch and Candelas jHj.) Thus, for £ = 1/6 we must keep subleading 

terms. This includes keeping the subdominant term in x, an d the distinction between 

t{z) and t = t(z = za/r), 

2r — a j r — a\ 2 _~ □ r — a , . 

X~t 1 , zt^zt-t 3 z , 22 

t r \ r J r 

as well as the next term in the uniform asymptotic expansion of the Bessel functions, 
s,(x) ~ -y/zte" 1 (1 + ui (t)/v + ...), ei(x) ~ Vrte~ vri (I -u x {t)/v + ...), (23a) 

4(x) ~ (1 + + ...), e[(x) ~ (1 - ^(t)/^ + ...), (235) 



where ui(f) = (3t - 5t 3 )/24, v^t) = (3t + 7t 3 )/24, as I -> oo. Including all this, it is 
straightforward to recover the well-known result (J2J of Deutsch and Candelas for strong 
coupling (Dirichlet BC). Following the same process for weak coupling, we find that the 
leading divergence in order n, 1 < n < 3, is (r a±) 

Ml "'-(l)" 1 44 1 0^ a (a-V" ( "- 1)( " + 2>r(3 -" ) - (M) 
Note that the subleading 0(A) term again vanishes. Both of these results apply for the 
conformal value £ = 1/6. 

The above results for the conformally coupled scalar show that the inverse linear 
divergences which occur in either order A or A 2 cancel between inside and outside, when 
one computes the total energy, while the divergence encountered at n = 3 is logarithmic: 

a-\ A 3 1 . , A 3 , r — a , . 

u { ) ~ -= o-r -> — n , 25 

a 7 144tt 2 V ; 1447T 2 a 7 a ' V ; 

where the latter form is shown by explicit calculation, rather than continuing in n. The 
integral of this, however, is finite, so this does not signal any difficulty. 
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5. Surface and shell energy 

However, as discussed first by Dowker, Kennedy and Critchley [T%| ITH] . and later 
elaborated by Saharian and Romeo [23 EE] , and put in a broader context by Fulling 
[22], for situations when other than Neumann or Dirichlet boundary conditions apply, 
an additional term must be supplied in calculating the energy, a term which resides 
entirely on the surface. For the case of the general stress tensor, that extra term is 



/ dS- VG(x,x') 



0i , - ■ (26) 

Z,l J x=x' 

where the direction of the normal is out of the region in question, which arises from the 
T°* component of the stress tensor. The total energy in a given region is not, therefore, 
just the integral of the local energy density, but has this extra contribution 



S = I (dr) (T 00 ) + £=jj (dr) J |^ 2 £(r, r^, (27) 



which is independent of £. (r is a point-splitting regulator [23].) The latter expression 
has a rather evident interpretation in terms of summing zero-point energies. The surface 
energy cancels for a nonsingular potential when computing the total energy in all space. 

In the limit of h — > oo for the region in the shell, a_ < r,r' < a + , the reduced 
Green's function becomes (for further details about this limit see 

1 ei{na)si{Ka) 



2nrr' 1 + -^ei(Ka)si{Ka) 



V Xh . .. \\h . . 
cosh (r — r ) + cosh (r + r — a + — a_ 



28 



In the thin shell limit (5 — > 0) this leads to an energy density in the shell nearly 
independent of r, leading to the energy (e = rg/a, y = \x\) 

Es = - 40 £(2/ + 1) r dx ^iffry, r e 1 - (29) 

However, we have to include the surface term (J2l)|) in the shell at r = a±, which exactly 
cancels this: £ s = E s + <B S = 0, because the total energy of the shell is given by (|27|) 
integrated over the volume of the shell, which clearly vanishes as the thickness of the 
shell 5 — > 0. However, we shall see shortly that E s plays a special role. 

6. Total energy of A sphere (5 = 0) 

Likewise, if one integrates the interior and exterior energy density, and includes the 
surface energy, one gets, for arbitrary £, the total energy as found by Bordag et al [13J, 



1 00 Z" 00 d 

£ = E in + E out + £ = --— £(2Z + 1) / dxy— In 

4na , n J-oo ay 



1=0 



1 + -I u (y)K v (y) 
a 



e ixe , (30) 



dy 

exactly that obtained from the integral (j2*7|) of the Green's function. 

However, there is more to say here. As noted above, the integral of the local 
energy, inside and outside the sphere, is finite perturbatively, because of cancellations 
between inside and outside, for the conformally coupled scalar. But it is well known 
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that divergences occur in the total energy at order A 3 . These evidently must arise from 
the surface term. So let us consider the latter, which is given in the outside region by 



00 1 r°° <\f 8 



e*™ |<| = «. (31) 

r=r'=a 



1=0 

In the strong coupling limit, there is, of course, no surface term. This is because then 



r,r' > a : gi(r, r' 



Sj(Kr<)ej(«r>) - Sl ha ei ( Kr )ei(nr') 



(32) 



nrr' 

which vanishes on the surface, and has a derivative proportional to the Wronskian. 

In general, in the thin-shell limit, the sum of the inside and outside surface terms 
is given by 

t= ^-^L^ + \!tSliv) ^ (33) 

Perhaps not remarkably, this is precisely the same as the integrated local shell energy 
E s (|29|). Thus the surface energies within and outside the shell regions cancel. (This is 
generally true, as follows from the continuity requirements on the Green's function.) 

For weak coupling, we expand this in powers of A. Perhaps the easiest way to 
isolate the asymptotic behavior is to use the leading uniform asymptotic expansion, 
I v (x)K v (x) ~ t/2u. This yields the following expression for the nth term in the total 
surface energy, (e = 0, analytically continued in n from Ken > 3) 

(-]) n ( a \ n r f —) 

ew ~-^ (i - 4 «y T^) {2 '" 2 - m "- 2) - <34) 

Note that this expression vanishes for n = 2; in this approximation the order A 2 term 
in the energy arises only from the local energy density. However, for n = 3 we obtain 
for the conformal value, £ = 1/6, 

«« - 5S>>. (35) 

precisely the divergent term in the energy first found by the heat kernel calculation of 
Bordag, Kirsten, and Vassilevich ^3]. The universality of these results supports the 
hypothesis of analyticity in the order n. Alternatively, if we keep e 7^ 0: 

(2) ^ 2 [°° dz I (3) A 3 . , nn . 

£ ( > ~ r / 7, = 0, <£ (3) ~ -hie. (36) 

24vrea 3 J -00 z z 2 + 1 ' 127ra 4 v ; 

The former integral vanishes by oddness, while the 0(A 3 ) term is logarithmically 
divergent as e — > 0. Thus, by expanding (jSOJ) in powers of A, £^ = A 2 /327ra 3 , is 
unambiguously finite, while £^ is unambiguously divergent. 



7. Conclusions 



For the case of a massless scalar field in a spherically symmetric step-function shell 
potential, we have shown that there is a net effective surface energy in the thin shell limit, 
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to be added to the integrated local energy density for the inside and outside regions, 
which is exactly the integrated local energy density of the shell. This shell energy, for 
the conformally coupled theory, is finite in second order in the coupling, but diverges in 
third order. We show that the latter precisely corresponds to the known divergence of 
the total energy in this order. Thus we have established the suspected correspondence 
between surface divergences and divergences in the total energy, which has nothing to do 
with divergences in the local energy density as the surface is approached. This precise 
correspondence should enable us to absorb such global divergences in a renormalization 
of the surface energy, and should lead to further advances of our understanding of 
quantum vacuum effects. Further details are given in [23]. 
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